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^ ' The recently published final results of experiment E04-108 in Jefi'erson Lab 

l/^ ' Hall C extend the recoil polarization measurements of the ratio of the pro- 

l/^ , ton electric and magnetic form factors to = 8.5 GeV^, an increase in 

00 ■ coverage of more than 50%. A global fit of and to selected data for 

C " 3 ' electron-proton elastic scattering cross sections and polarization observables is 

' presented, illustrating the statistical impact of the new results. 
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The elastic electromagnetic form factors of the nucleon are among the most 
fundamental quantities describing its structure. In the one-photon-exchange 
or Born approximation, the helicity-conserving Dirac form factor FiiQ"^) 
and the helicity-flip Pauli form factor F2{Q^) fully characterize the ampli- 
tude for elastic electron-nucleon scattering at a spacelike momentum trans- 
fer = —q^ > 0. These form factors are accessible experimentally through 
both cross section and double-polarization measurements, and are often ex- 
pressed in terms of the equivalent linear combinations Ge ^ Fi— TF2 (elec- 
tric) and Gm — Fi+ F2 (magnetic), which are more readily extracted from 
experimental observables. The Born differential cross section for elastic ep 
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scattering is given by^ 

d^e Q^\Ej 2 1 + r 

where a is the fine-structure constant of electromagnetism, E,, and E'^ 

arc the lab incident and scattered electron energies, 6^ is the lab elec- 
tron scattering angle, r = Q^/AM"^, M is the nucleon mass, and e = 
[l + 2(1 + r) tan^ ^] corresponds to the longitudinal polarization of the 
virtual photon in the single-photon-exchange picture. 

The polarization of protons scattered elastically by longitudinally po- 
larized electrons has nonvanishing longitudinal and transverse components 
relative to the momentum transfer in the scattering plane given by^'^ 

loPi = V^7(TT7) tan^ I^L+^Gi,, (2) 

loPt = -2VT(1 + T)tan jGeGm, (3) 

where /q = G% + ■jG\[. The component P„ normal to the scattering plane 
is exactly zero in the Born approximation. The ratio Pt/Pi is directly pro- 
portional to Ge/Gm'- 

Pt Ee+E'^ ^ e. 

The ratio (4) provides enhanced sensitivity to at large Q^, where Gm 
dominates the cross section (see equation (1)). Precise recoil polarization 
experiments'*"^ first established that R decreases rapidly with > 1 GcV^, 
in strong disagreement with cross section measurements. The extraction of 
a small G\ term compared to a dominant jG\j term from L/T separation 
measurements becomes highly sensitive to incompletely understood higher- 
order QED effects at large Q^, including two-photon exchange (TPEX).^ 
The enhanced sensitivity of the recoil polarization method to G^; and the 
typically smaller relative importance of radiative corrections and TPEX to 
the ratio (4) have led to a general consensus that the polarization data most 
reliably determine G^ at large Q^. 

Extending the accurate measurements of the proton and neutron elec- 
tric and magnetic form factors to the highest accessible provides cru- 
cial experimental input to understanding the transition between the non- 
perturbative and perturbative regimes of QCD. In addition to constraining 
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QCD-inspired phcnomcnological models, precise form factor data provide 
important model-independent information about the structure of the nu- 
cleon. The recently published results of experiment E04-108^ in Jefferson 
Lab's Hall C extended the recoil polarization data for RtoQ^= 8.5 GeV^. 
In this paper wc illustrate the statistical impact of the new data in a global 
fit to elastic electron-proton scattering data, including both cross section 
and polarization measurements, using the Kelly parametrization^ of 
and G^^. 

The Kelly parametrization of the miclcon form factors^ is given by the 
ratio of a polynomial of order n to a polynomial of order n + 2 in r = 
gV4M2: 

This parametrization satisfies G(Q^) cx as ^> oo and gives the ap- 
propriate static (Q^ 0) limit if one fixes oq = 1 for G^ and G^j/ jjLp. This 
parametrization is chosen to have the large- behavior expected from di- 
mensional scaling laws in perturbative QCD.^"'^^ The original 2004 analysis 
by Kelly using this parametrization found that a four-parameter fit with 
n = 1 could describe the data for G^ and G^" with a reasonable x^. For 
simplicity's sake, the results presented here use the same choice. 

The published elastic ep scattering data chosen for the fit include 421 
differential cross section measurements spanning 0.005 GeV^ < < 
31.2 GeV^ and 53 polarization measurements spanning 0.162 GeV^ < < 
8.49 GeV^. The differential cross section data were taken from Refs. 12-23, 
while the polarization data were taken from Rcfs. 5,6,8,24-30. While by no 
means comprehensive, the chosen data are sufficiently representative of the 
precision and coverage of all elastic ep scattering data. All data included 
in this analysis were obtained from hydrogen targets. 

The best fit parameters for G^ and G^ were obtained in a simultaneous 
global analysis of cross section and polarization data by minimizing the 
function defined as 



'-^"R,data 



E 



^p,data 



(6) 



where Rp — /ipG^/G^, G^ and G^/iJ,p are given by the formula (5) with 
n = 1 and oq = 1, of and bf'^ are eight adjustable parameters to be 

determined, (TR^data = e(^ + T)adata/ cfmou (sec equation (1)), and AaR^data 
and ARp^data are the uncertainties in the experimental data. For the cross 
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section data, statistical uncertainties and overall normalization uncertain- 
ties, when quoted with the published results, were included in AaR^data- 
For the polarization data, only statistical uncertainties were included. 
The inconsistency between the Rosenbluth and polarization data for 
> 1 GcV^ was handled using an iterative procedure. Three iterations 
of the fit were performed. For all cross section data with > 1 GeV^, the 
value of in equation (6) was replaced by a fixed value of calculated 
from the results of the previous fit, or, in the case of the first iteration, using 
the results of Kelly's 2004 analysis.^ This replacement has the effect that 
for >1 GeV^ , is entirely determined by polarization data, while G^ 
is determined by cross section data, with fixed by the polarization data 
for Rp. For < 1 GeV^, cross section and polarization data are treated 
on an equal footing for both G^ and G^. The numerical minimization was 
carried out using the TMinuit class within the ROOT libraries. Typically, 
no significant improvement of the fit was found after two iterations of the 
starting function for G^. 

This fit procedure was performed for two data sets identical in every 
respect except that the new JLab Hall C data for Rp were excluded from 
the first data set and included in the second set. Table 1 shows the results 

Table 1. Fit results for and G^/fip, with ("Old") and without ("New") 
the data of Rcf. 8. Goo/Gd = 16Afp/A^ (ai/fe) is the asymptotic value of 
C^/Go and G'^ / {fj,pG o) ■ The "standard" dipole form factor is defined as 
Gd = (1 + Q2/a2)-2^ ^ith a2 = 0.71 GeV2. 





Old G| 


Old gIj 


New 


New Glj 


Ol 


-.390 ± .081 


.075 ± .023 


-.299 ± .072 


.081 ± .023 


bi 


11.01 ± .10 


11.14 ± .08 


11.11 ± .09 


11.15 ± .08 


b2 


13.57 ± .84 


18.44 ± .16 


14.11 ± .83 


18.45 ± .16 


63 


11.4 ± 4.5 


5.12 ± .65 


15.7 ± 4.2 


5.31 ± .66 


Goo /Go 


-.84 ± .49 


.360 ± .064 


-.47 ± .23 


.376 ± .061 


QHGl=0) 


9.0 




11.8 





of the fit with and without the new Jefferson Lab results. The parameter 
errors quoted are the standard errors calculated from the diagonal elements 
of the c;ovariance matrix. Significant correlations among the parameters for 
each form factor were found. These correlations imply that the uncertain- 
ties in the asymptotic form factor values and the uncertainties in G^ and 
G\[ at any given cannot be calculated directly from the parameter er- 
rors, but instead require the full covariancc matrix. On the other hand, the 
correlation coefficients between the parameters describing G^ and those 
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describing were generally small, at the level of a few percent. Table 
2 shows the asymptotic values of Q^i^f ^^'^ Q'^F^ for the fits with and 
without the new high data. 



Table 2. Asymptotic values of Q*Fi and Q^F2, with and without 
Ref. 8. 



Old Q^Fl^ 


New Q^Fl^ 


Old Q*Fl^ New Q'^Fl^ 


3.28 ± .93 


2.69 ± .51 


.507 ± .090 .529 ± .086 



Table 3 illustrates the quality of the fit in terms of x^- Although the 
overall of the fit x^/n.d.f. = 828.4/466 = 1.78 might be regarded as 
relatively poor, the extent to which the reduced x^ exceeds one per datum 
reflects the inconsistency between the Rosenbluth and polarization meth- 
ods, as shown in table 3. The polarization data and the cross section data 
for <1 GeV^ are consistent with the fit results for and G^^ . Relative 
to the number of data points, the cross section data for 1 < < 10 GeV^ 
contribute disproportionately to x^ ) since the e-dependence of the measured 
cross sections disagrees with the slope predicted by the polarization data. 



Table 3. of the final global fit with new 
data included. 



Observable 






o-fl, < 1 GcV^ 


269 


293.8 


o-R, 1 < < 10 GcV^ 


142 


483.3 


cTfl, Q2 > 10 GeV^ 


10 


3.96 


allQ2 


53 


47.4 


All data 


474 


828.4 



Figures 1 and 2 show the fit results together with the experimental 
data. In figure 1(a), polarization data for Rp were converted to G^ data 
using the fit result for G^ and appropriate error propagation. The best fit 

curves are shown with standard Icr uncertainty bands calculated from the 
full covariance matrix. As shown in figure 1(a), the new G^ data* favor 
a slowing rate of decrease of G^ with compared to previous data, and 
shrink the uncertainty in G^- for > 5 GeV^ by roughly a factor of two. 
The improved constraint on in the high-Q^ region shown in figure 1(b) 
is similarly dramatic. 

Figure 2 shows the results for G^^/npGo and Q'^F^. The interesting 
feature of the G^^ results is that the fit is systematically higher than the 
published data in the region of the discrepancy. Since the cross section 
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(a) GIIGd (b) Q^Fl 

Fig. 1. Fit results for G^^ j G oia) , and Q^Fj (b), before (red dot-dashed) and after (blue 
dashed) including the results of Ref. 8, with Icr uncertainty bands. Published values of 
extracted from cross section (open circles) and polarization (filled circles) measurements. 
Results of Ref. 8 (magenta squares). 

measurements either extracted or assumed i?p « 1, the contribution of G\ 
to (JR was overestimated, lowering the extracted values of G\j . This result 
for is quantitatively and qualitatively similar to the result of a recent 
global analysis including TPEX corrections to the cross section data.'^^ In 
contrast to this analysis, which neglects TPEX completely, the authors of 
Ref. [32] used the Kelly parametrization with n = 3 and achieved a sub- 
stantially better x^, since the applied TPEX corrections largely reconciled 
the discrepancy. Although the new Rp data barely affect G^^ and its un- 
certainty, the effect on F\ is more significant, an entirely expected result 
given the dramatic improvement in and the definition Gm = F\ ^ F^. 
For yet larger Q^, the fit result for G\^ again agrees with the published 
data, reflecting the dominance of cta by G\^^ shown in figure 3. 

In conclusion, new high measurements of R — fj,pG^/G^ in Jefferson 
Lab's Hall C have significantly extended the range of for which G^ is ac- 
curately determined. The impact of these new data on the world database of 
proton form factors was studied in the context of an empirical parametriza- 
tion with physically reasonable assumptions for low and high asymp- 
totic behavior; i.e., G^(0) = 1 and G^j(O) ~ /ip, and ^f(2) Q^^^^-* as 
oo. In this context, the new data for Rp reduce the uncertainty in 
G^ by roughly a factor of two in the measured region and in the ex- 
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Q^ GeV^ , GeV' 



(a) Gl,/i^^r,GD) (b) Q^Ff 

Fig. 2. Fit results for G\,j / {^ipG £i){a) and Q'^F^{h), before (red dot-dashed) and after 
(blue dashed) including the results of Ref. 8. Published values of G^^ extracted from 
cross section (open circles) measurements. 




Fig. 3. G|,/(G|, + rG\j), calculated from the results of the global fit, represents the 
maximum fraction of the reduced cross section carried by Ge ■ 

trapolation to higher Q^. Looking to the future, planned recoil polarization 
measurements of i?p to ~ 15 GeV^ and precision elastic ep differential 
cross section measurements to ~ 18 GeV^ using JLab's upgraded 11 
GeV electron beam will complete the experimental knowledge of proton 
electromagnetic form factors in the spacelike region attainable with present 
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day accelerators. 
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